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Abstract 

We report the first observation of Z/ 7 * production in Compton scattering of quasi-real 
photons. This is a subprocess of the reaction e’''e“ —> e+e“Z/ 7 *, where one of the final state 
electrons is undetected. Approximately 55 pb“^ of data collected in the year 1997 at an e'''e“ 
centre-of-mass energy of 183 GeV with the OPAL detector at LEP have been analysed. The 
Z/ 7 * from Compton scattering has been detected in the hadronic decay channel. Within well 
dehned kinematic bounds, we measure the product of cross-section and Z/ 7 * branching ratio to 
hadrons to be (0.9±0.3±0.1) pb for events with a hadronic mass larger than 60 GeV, dominated 
by (e)eZ production. In the hadronic mass region between 5 GeV and 60 GeV, dominated by 
(e)e 7 * production, this product is found to be (4.1±1.6±0.6) pb. Our results agree with the 
predictions of two Monte Carlo event generators, grc4f and PYTHIA. 
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1 Introduction 


Single Z boson production in e+e“ collisions [Q], e’''e“ —> (e)eZ, will be the dominant source of 
Z bosons at linear e’''e“ colliders with centre-of-mass energies above 500 GeV P|. At LEP2, 
however, the cross-section is nearly two orders of magnitude below the “radiative return” process 
e’''e“ —yZ. The elementary subprocess of single Z boson production isQ e “7 e~Z where 
a quasi-real photon radiated from one of the beam electrons scatters off the other electron 
producing a Z as shown in Fig. [^. This process is ordinary Gompton scattering with the outgoing 
real photon replaced by a Z. Its cross-section has hrst been calculated for ep collisions 0, where 
the incoming photon is radiated from the proton. Another variant of Gompton scattering is 
e “7 —> e“ 7 *, where this time the outgoing real photon is replaced by a virtual one. The 
observable hnal state (e)eff will in both cases be the scattered electron e and a fermion pair ff 
from the Z/y* decay, while the other electron (e) usually remains unobserved in the beam pipe. 



Figure 1: Diagrams for the process e^e —(e)eZ/ 7 *. 
^Charge conjugation is implied throughout the paper except when otherwise stated. 
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The final state (e)eff may contribute an important background to many searches for new 
physics, for example final states involving large missing masses (ff = uu), leptoquark production 
in deep inelastic electron-photon scattering (ff = qq), new particles decaying to Z bosons, and 
excited electrons in the decay channel eff. Excited electrons can for example directly contribute 
in the s and u channel diagrams of the signal in Fig. |^. 


For Z boson production in Compton scattering e{k)'y{p) —>• e{k')Z{p') of real photons 
= 0), the cross-section dependence on the Mandelstam variables s = {k'+p')^, t = {k' — kY, 
u = {p' — kY is P 

da 1 (u 2Myi s\ . . 

dt s-^ \s us uj 

For Mz = 0 the well-known terms for ordinary Compton scattering remain. The typical trans¬ 
verse momentum scale of the scattered Z/y* bosons is small, leading to a singularity at h = 0 
of the virtual electron propagator in Fig. ^b. For incoming quasi-real photons (p^ ~ 0) in ep or 
e’''e“ collisions, the dominant regulating effect for this divergence is not the electron mass, but 
small, non-zero, incoming photon masses squared p^, which introduce via the replacement 0 


u 


u + p^ — 
s 


( 2 ) 


a lower cutoff in the denominator of Eq. (|^). A simple equivalent photon approximation (EPA), 
as used in Ref. [Q, where the integration over the range of small photon virtualities leads to 
an effective on-shell incoming photon flux, overestimates the cross-section by about a factor of 
two [^. In order to properly describe the process, the p^ spectrum of the incoming photons has 
to be retained fully, or a modified EPA has to be introduced. In any case the results will 
be sensitive to the modelling of the p^ spectrum. In this paper the theoretical expectations are 
represented by Monte Carlo event generators which use different approaches for obtaining the 
p^ spectrum of the incoming photons. These are compared with experimental data for the first 
time. 


We have analysed the process e+e“ —>• (e)eZ/ 7 * using data collected with the OPAL detec¬ 
tor at an e'''e“ centre-of-mass energy of about 183 GeV, corresponding to a total integrated 
luminosity of approximately 55 pb~^. The Z or the 7 * is observed via its decay to hadrons. 


2 The OPAL detector 


A detailed description of the OPAL detector can be found elsewhere [|]. Subdetectors which 
are particularly relevant for the present analysis are described here briefly. The central detector 
consists of a system of tracking chambers providing charged particle tracking over 96% of the full 
solid angle inside a 0.435 T uniform magnetic field parallel to the beam axi^. Starting with the 
innermost components, it consists of a high precision silicon microvertex detector, a precision 
vertex drift chamber, a large volume jet chamber with 159 layers of axial anode wires and a set 
of 2 ; chambers measuring the track coordinates along the beam direction. With the jet chamber, 

^The OPAL coordinate system is a right-handed system which is defined such that the z-axis is in the 
direction of the electron beam and the a;-axis is horizontal and points towards the centre of the LEP ring; 6 is 
the polar angle with respect to 2 and (j) is the azimuthal angle about the z axis. 
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the momenta of tracks can be measured with an accuracy of cfp/p^ = 2 x 10 ^ GeV ^. The jet 
chamber also provides energy loss measurements which are used for particle identihcation. 

A lead-glass electromagnetic calorimeter (ECAL) located outside the magnet coil covers 
the full azimuthal range with excellent hermeticity in the polar angle range of | cos 6^1 < 0.82 
for the barrel region and 0.81 < | cos 6^1 < 0.984 for the endcap region. The forward detectors 
and silicon tungsten calorimeters 0 located at both sides of the interaction point measure 
the integrated luminosity and complete the geometrical acceptance down to 33 mrad in polar 
angle. The magnet return yoke is instrumented with streamer tubes for hadron calorimetry 
and is surrounded by several layers of muon chambers. 


3 Event simulation 


The e’''e“ ^ (ejeZ/y* signal events are generated using the grc4f Monte Carlo program |^. 
Hadronization is performed by JETSET |^. A sample of approximately 170 times the integrated 
luminosity of the data was used. In grc4f, an automatic computation system for Feynman 
diagrams, GRACE is invoked, yielding the differential cross-sections of tree-level diagrams 
by performing a calculation of the matrix elements and a Monte Carlo integration over the phase 
space. Initial state radiation corrections are also implemented. The delicate distribution of 
the incoming photons for the ey Compton scattering is thus obtained by a numerical integration 
of the matrix element. 


As a second Monte Carlo event generator for the process e'''e —>■ (e)eZ/y* we used the 
program. To include the full phase space for this low-pt process, we lowered 


PYTHIA M 


the pt cutoff in the ey rest frame from the default of 1 GeV to 10 MeV in the generation 
of the events, as suggested in the PYTHIA manual. Without this change PYTHIA would 
underestimate largely the cross-section. In PYTHIA, the signal process is simulated according 
to Eqs. (|I|) and (|^). An equivalent flow of incoming photons is modelled with a p^ distribution 
obtained from a hrst order QED expression [|^]. This approximation implemented in PYTHIA 
is used as a cross-check, with a sample size of about 1/10 of that generated with grc4f. 

Backgrounds from other four-fermion final states have been studied using several samples 
also generated with grc4f and JETSET. Background contributions from two-photon processes, 

for events with low 


e“'"e 


e' e 


hadrons, were evaluated using PHOJET [P 


and 

HERWIG [O for high-Q^ processes. Double counting of the eeqq hnal state between signal. 


two-photon and other four-fermion samples is avoided by separating the contributions from 
different classes of diagrams, or by performing kinematic cuts. Two-photon production of 
was modelled with the VERMASEREN |T^ generator. Other background processes involving 
two fermions in the hnal state were simulated using PYTHIA for the channel e’''e“ —> qq and 
KORALZ 


141 for the channel e'''e 


T^T . Background contributions from other processes 


have been checked and found to be negligible. 


All Monte Carlo events were passed through the full simulation of the OPAL detector [p 
and then subjected to the same reconstruction and analysis procedures as the real data. 
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4 Signal definition 


In most of the e’''e“ —> (e)eZ/ 7 * events the hrst incoming electron, which radiates the photon 
taking part in the Compton process, remains undetected in the beam pipe due to the low mo¬ 
mentum transfer In contrast, the second electron is often emitted at large angles. Therefore 
the predominant signature of the process is a single scattered electron isolated from a hadronic 
system containing one or two hadron jets from the Z/ 7 * decay. The differential distribution 
is peaked at low (low mass qq systems) and at low \u\ where the outgoing Z/ 7 * is emitted 
almost in the direction of the second incoming electron. Since for low p'^ and low |m| a huge 
part of the cross-section remains unobserved in the beam pipe, it is not practical to measure 
the signal process e’''e“ —>■ (e)eZ/ 7 * in the whole phase space corresponding to the diagrams 
in Fig. 2. Furthermore, considerable uncertainties due to hadronic resonances for low mass 
7 * provide a supplementary reason for limiting the phase space. The e’''e“ —> (e)eZ/ 7 * signal 
process is therefore dehned by the diagrams in Fig. ^and by applying the following additional 
kinematic cuts on the eeqq hnal state. 


• One of the primary electrons is required to remain in the beam pipe: | cos 6 *eJ > 0.9995. 
This angular cut is chosen close to the edge of the silicon tungsten calorimeter. 

• The second beam electron has to form a large angle with respect to the beam: | cos 9^^ \ < 
0.985, corresponding to the edge of the electromagnetic endcap calorimeter. 

• In order to detect the hadronic hnal state from the Z/ 7 * decay, we demand at least one 
primary quark with | cos^ql < 0.985. 

• We require the invariant mass of the primary quark pair to be greater than 5 GeV. 


The cross-section times branching ratio predicted for e’''e“ —> (e)eZ/ 7 * events with an eeqq 
hnal state satisfying the signal dehnition has been evaluated using the grc4f and PYTHIA 
programs at an e’''e“ centre-of-mass energy of 183 GeV. We obtain values of (4.80±0.02) pb 
for grc4f and (4.76T0.04) pb for PYTHIA. 

In addition to the ey Gompton scattering in the process e’''e“ —> (e)eZ/ 7 *, there are two 
other classes of processes which lead to an eeqq hnal state, namely e’''e“ annihilation, dominated 
by Z/ 7 * — Z/ 7 * intermediate states, and two-photon interactions. Of the former process, only 
a negligible cross-section of 0.002 pb fulhls the kinematic cuts of our signal dehnition, since the 
topology with exactly one electron very close to the beam direction and the other electron visible 
at large angles is not preferred. Among the two-photon interactions, on the other hand, exactly 
this “single-tagged” topology is typical for deep inelastic ey scattering with large momentum 
transfer squared]^, = —Q^. According to PHOJET and HERWIG, two-photon events which 
satisfy the above kinematic cuts are produced with a cross-section of approximately 5.6 pb. 

Using grc4f we checked that possible interferences between the diherent classes of Feynman 
diagrams contributing to the eeqq hnal state, ey Gompton scattering, e’''e“ annihilation, and 
two-photon interaction, are negligible. It is therefore justihed to treat each process separately 

^Note, that in two-photon processes is by definition identical to t in our signal process. 
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and to add the predictions of the respective Monte Carlo generators. Only e’''e“ —(e)eZ/ 7 * 
events which fnlhl the signal dehnition are treated as “signal”. All other processes leading 
to an eeqq hnal state are treated as “backgronnd” even if the kinematic signal dehnition is 
fnlhlled. Likewise, e’''e“ —*• (e)eZ/ 7 * events that do not fnlhl the signal dehnition are treated 
as “backgronnd”. 


5 Event preselection 


We perform a preselection of events with a hadronic hnal state and an isolated electron can¬ 
didate. Only tracks and clnsters which satisfy standard qnality criteria are considered in the 
analysis. The following preselection cnts are applied. 


• We reject low mnltiplicity events by reqniring the snm V^uit of the nnmber of tracks and 
electromagnetic clnsters not associated with any track to be greater than 8. 


• We search for electron candidates among all tracks nsing the following criteria. We reqnire 
the probability of the measnred ionization energy loss dE /dx of the track being consistent 
with an electron hypothesis to be greater than 5%. Fnrthermore we reqnire Eg/pe > 0.6 
where Eg is the total energy in the electromagnetic calorimeter associated with the track 
and Pe is the track momentnm. For all electron candidate tracks we dehne a cone energy 
Egone as the snm of all track momenta and electromagnetic clnster energies in a cone with 
a half angle of 0.2 rad aronnd the track, exclnding the momentnm and the clnster energy 
associated to the electron candidate track. If more than one candidate is fonnd in an 
event, the one which has the smallest i^cone is chosen as the signal electron. To snppress 
events with fake electron candidates, we reqnire the angle between the signal electron and 
the closest track to be greater than 0.25 rad and Eg to be larger than 2 GeV. 


In order to be consistent with the signal dehnition, we reqnire the invariant mass of the 
qq system rriqq to be larger than 5 GeV. This mass is calcnlated from the momentnm 
of all tracks and the energy deposited in the electromagnetic and hadronic calorimeters, 
the forward detectors and silicon tnngsten calorimeters, exclnding those associated with 
the signal electron, nsing an algorithm which corrects for donble connting of energy 
between tracks and calorimeter clnsters. 


In Fig. ^ we compare the measnred distribntions of mqq and Eg after preselection with 
the corresponding Monte Garlo expectations. The gross featnres of the data are in satisfactory 
agreement with the prediction at this early stage of the analysis. Discrepancies, especially at low 
valnes of mqq and Eg, can be attribnted to systematic nncertainties in the two-photon process 
e'''e“ —> e’''e“ -I- hadrons which is difhcnlt to model in this kinematic region. The apparent 
excess of Monte Garlo events in the range 105 GeV < mqq <160 GeV has no inflnence on the 
analysis, since it is populated nearly exclusively by background events, which will be rejected 
later on. A total of 1558±25 events is expected, where the main contributions are from two- 
photon interactions, e'''e“ —> and multihadronic events. In the data 1882 events are 

found after the preselection (see Table j^). As seen from Fig. H, the data excess is located in 
kinematic regions where the two-photon background dominates. 
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The expected signal contribution is 38.0T0.5 e’''e“ —>• (e)eZ/ 7 * events, corresponding to an 
average signal efficiency of (14.5±0.2)% obtained with the grc4f Monte Carlo program. The 
corresponding value obtained from PYTHIA is (16.7±0.6)%. A large part of the preselection 
inefficiency for all Z/ 7 * masses is due to the typical small transverse momenta of the electron 
and the Z/ 7 * in the hnal state. Signal electrons which form large angles with respect to the 
beam therefore often have momenta below 2 GeV. In addition, for low 7 * masses, the 7 * decay 
products tend to be in a rather narrow cone preferring forward directions, where particles may 
escape undetected along the beam pipe. As a result many of these events fail the multiplicity 
cut. This reflects in a signihcantly lower preselection efficiency for smaller hadronic masses. 
Considering only signal events with rriqq < 60 GeV, grc4f predicts the preselection efficiency to 
be (9.5±0.2)%, while for m^q > 60 GeV the efficiency is (50.3±0.7)%. The PYTHIA prediction 
is consistent with these numbers. 


Cut 

Number of expected events from MC 
(e)eZ/7* 77 qq 4 f other Sum 

OPAL 

data 

Preselection 

38 . 0 ± 0.5 

1251 . 9 ± 23.9 

113 . 8 ± 3.6 

143 . 1 ± 0.9 

11 . 6 ± 0.6 

1558 . 4 ± 24.7 

1882 

— Qe • COS0qq > 0.75 

33 . 0 ± 0.4 

575 . 0 ± 16.5 

42 . 6 ± 2.2 

11 . 6 ± 0.3 

3 . 5 ± 0.4 

665 . 7 ± 16.7 

766 

—Qe • cos A < 0.75 

26 . 5 ± 0.4 

209 . 5 ± 10.1 

30 . 8 ± 1.9 

9 . 3 ± 0.2 

2 . 0 ± 0.3 

278 . 1 ± 10.3 

294 

1 COS ^miss 1 ^ 0.95 

24 . 3 ± 0.4 

128 . 9 ± 7.8 

20 . 2 ± 1.5 

2 . 9 ± 0.1 

1 . 1 ± 0.2 

177 . 5 ± 7.9 

193 

Pmiss ^ 30 GgV 

21 . 6 ± 0.3 

21 . 0 ± 2.0 

15 . 9 ± 1.4 

2 . 0 ± 0.1 

0 . 7 ± 0.1 

61 . 2 ± 2.5 

74 

electron isolation 

17 . 5 ± 0.3 

8 . 2 ± 1.2 

1 . 9 ± 0.5 

1 . 4 ± 0.1 

0 . 6 ± 0.1 

29 . 6 ± 1.3 

36 

< 35 GeV 

16 . 7 ± 0.3 

3 . 1 ± 0.8 

1 . 8 ± 0.5 

1 . 4 ± 0.1 

0 . 3 ± 0.1 

23 . 2 ± 0.9 

27 


Table 1: Numbers of expected and observed events for a luminosity of 55 pb“^ after application 
of each cut. Only the most important sources of background are shown separately. The number 
of expected signal events is obtained using the grc4f generator. The numbers of background 
events are evaluated using the Monte Carlo samples (MC) described in the text. The errors are 
statistical only. 


6 Selection of e+e —(e)eZ/ 7 * events 


In order to reduce the remaining background six further cuts have been applied which can be 
divided into three groups. The hrst group consists the following angular cuts. 

• In the e’''e“ ^ (e)eZ/ 7 * signal process, the Z/ 7 * is usually scattered very close to the 
beam direction of the detected signal electron, given by —Qe • cos 6 * = 1 , where Qe is the 
charge (in units of e) of the signal electron. We therefore require —Qe • cos 6 *qq > 0.75, 
where cos^qq is the cosine of the polar angle of the qq system (see Fig. ^). 

• To suppress the otherwise irreducible background of single-tagged two-photon events, 
which fulhl the kinematic signal dehnition, we take advantage of the fact that the distri¬ 
bution of the polar angle cos 6e of the tagged electrons peaks around their beam direction. 
We thus only retain events where —Qe • cos^e is less than 0.75. This cut is shown in 
Fig. ^b. Note that two-photon events with fake electrons, distributed symmetrically in 
—Qe • cos 6 *e, largely dominate this spectrum. Gorrectly identihed tagged electrons of high 

two-photon events are visible as an asymmetry in the region near -|- 1 . 



















Signal events are farther characterised by an nndetected electron along the beam direction, so 
that cuts on the direction and amount of the missing momentum in the event can be applied. 


• We require the absolute value of the cosine of the polar angle of the missing momentum 
I cos^missl to be larger than 0.95 (see Fig. ^jc). The largest relative suppression of this cut 
occurs for e^e“ —>■ W'''W“ events, where one W boson decays to hadrons and the other 
to an electron and a neutrino, since the neutrino is not preferentially emitted close to the 
beam direction. 

• The distribution of the missing momentum after the cut on | cos is shown in Fig. 
Requiring a missing momentum of at least 30 GeV largely suppresses remaining back¬ 
ground with fake or conversion electrons from untagged two-photon events where both 
beam electrons stay under small angles with respect to the beam axis and tend to balance 
the momentum. 


The third group of cuts is applied to suppress events with fake electron candidates. 


• The minimum angle between the electron and the nearest charged track was required to 
be 0.65 rad. This cut removes most of the remaining multihadronic background as shown 

in Fig. 1^. 

• The sample still includes a small amount of two-photon background tagged by an electron 
in the forward calorimeters where the signal electron is fake or due to a double-tagged 
event topology. It is suppressed signihcantly by requiring the total energy deposit in 
the forward detectors and silicon tungsten calorimeters, F^fwd, not to exceed 35 GeV (see 
Fig. |f). 


The cuts, together with the expected and observed numbers of events, are shown in Table [^. 
Already after the second cut there remains no signihcant excess in the data compared with the 
expectation. The six selection cuts improve the signal to background ratio by about a factor 
of 100 compared with the preselection, while retaining about half of the signal events. 

The efficiencies after all cuts are given in Table Gonsidering only signal events with 
mqq < 60 GeV, dominated by (e)e 7 *, grc4f predicts an overall selection efficiency of (4.1±0.1)%. 
For TTiqq > 60 GcV, dominated by (e)eZ, the overall efficiency is (22.3±0.6)%. The correspond¬ 
ing numbers from PYTHIA, (4.3±0.3)% and (22.3±1.7)%, are consistent with grc4f within their 
statistical errors. From the two lowest lines in Table |^, however, it is evident that PYTHIA 
predicts a slightly different mix of high mass “eeZ”-type and low mass “ee 7 *”-type events than 
grc4f. This causes about a 15% relative difference between the two Monte Garlo generators for 
efficiencies averaged over the total ruqq range and explains the differences in the total preselec¬ 
tion efficiencies given at the end of Section Since the predicted efficiency of each single cut 
is consistent between the two simulations, this difference remains essentially unchanged after 
the hnal selection cuts. 

For calculating cross-sections we avoid assumptions on the mix of “eeZ”-type, and “eey*”- 
type events. We instead calculate the cross-sections separately, for masses above and below 
60 GeV, the point where the differential mqq distribution has its minimum, using the efficiencies 
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obtained from grc4f. After correcting for feed-throngh from the low to the high mass region, 
and vice versa, we obtain the resnlts listed in Table The measnred cross-sections agree with 
the predictions of the two Monte Carlo generators. 



5 GeV < mqq < 60 GeV 
“eey*” 

^qq > 60 GeV 
“eeZ” 

efficiency (%) 

4.1±0.1 

22.3T0.6 

signal expected 

9.5±0.2 

7.2T0.2 

background 

4.4±0.8 

2.1T0.4 

feed-through from 
neighbouring mass region 

0.2 

0.1 

OPAL data 

14 

13 

cross-section (measured) (pb) 

4.1±1.6 

0.9±0.3 

cross-section (grc4f) (pb) 

4.21±0.02 

0.59±0.01 

cross-section (PYTHIA) (pb) 

4.04±0.03 

0.72±0.02 


Table 2: Comparison of the OPAL data with the cross-sections predicted by the grcff and 
the PYTHIA Monte Carlo generators for 5 GeV < mqq < 60 GeV (“ee'^*” region) and for 
> 60 GeV (“eeZ” region). The efficiencies have been obtained with grc4f. Only statistical 
errors are given. The errors on the feed-through are negligible. 

In Fig. § we have plotted after all cnts the same distribntions as in Fig. The respective 
contribntions of the 7 * and of the Z to the measnred process are clearly visible in the mqq 
distribntion. The distribntion shows the preference for small electron momenta dne to the 
low transverse momenta involved in the ey scattering. 


7 Systematic studies 


Several sonrces of systematic errors have been considered. Uncertainties in the description of 
the signal process may canse a systematic error on the efficiencies. These come mainly from 
imperfect modelling of the detector and have been analysed by comparing the Monte Garlo 
simnlation to real data for the process e+e“ —> W’''W“ where one W decays hadronically and 
the other one to an electron and a nentrino. This process has the same observable hnal state 
as e+e“ —>■ (e)eZ/ 7 * and can therefore be used to check for any discrepancies in the description 
of the cut variables which would lead to an error in the selection efficiency. After selecting 
the e’''e“ —>• ^ qqe^Ue events using the procedure described in Ref. |1T7| . 


each cut 


was applied separately to this sample. The difference in selection efficiency between data and 
the Monte Garlo simulation was taken as the systematic error after subtracting quadratically 
the statistical error of the WW sample. In those cases where the statistical error was larger 
than the efficiency differences, we have conservatively taken the statistical error as systematic 
uncertainty. The results of the check are shown in Table In our error estimate we have 
considered the fact that some distributions of the process e+e“ —> ^ qqe“z/e are 

signihcantly different from those of the signal (e.g. —Qe-cos 6 'qq or | cos 6 *iniss I) • These systematic 
uncertainties have been used for both the “ee 7 *”-like and the “eeZ”’-hke hadronic mass region, 
since the efficiencies of each selection cut are similar in the two regions. As mentioned above. 
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the relatively smaller overall efficiency for “ee 7 *”-like events is due to the multiplicity cut in the 
preselection. To assess the systematic uncertainty of this cut, we have calculated the change in 
efficiency when the simulated A^muit distribution is shifted by ±1. Half of the resulting change 
is quoted as uncertainty for this cut in the hrst row of Table ^ The total relative systematic 
uncertainty on the signal efficiencies due to event simulation is then found to be 0.086 for high 
mass and 0.107 for low mass hadronic systems, respectively. 



5 GeV < mqq < 60 GeV 

> 60 GeV 


“eey*” 

“eeZ” 

^mult > 8 

0.063 

0.000 

— Qe • COS^qq > 0.75 

0.018 

0.018 

— Qe • COS^e < 0.75 

0.037 

0.037 

COs6*miss > 0.95 

0.047 

0.047 

Pmiss > 30 GeV 

0.036 

0.036 

electron isolation 

0.046 

0.046 

h/fwd < 35 GeV 

0.009 

0.009 

event simulation 

0.107 

0.086 

generator distributions 

0.063 

0.063 

background 

0.070 

0.018 

Total 

0.142 

0.108 


Table 3: Relative systematic uncertainties of the cross-section measurements. The entry “event 
simulation” is the quadratic sum of the signal efficiency uncertainties for the single cuts listed 
in the rows above it. 

A smaller contribution is the uncertainty of the efficiency due to systematic errors in the 
Monte Carlo generator distributions. This has been estimated from the relative difference 
between the measured “eey*” and “eeZ” cross-sections obtained using the efficiencies from 
grc4f and those from PYTHIA. The relative cross-section differences are very similar in both 
hadronic mass regions, and statistically compatible with zero. Adding the cross-sections of both 
regions, a relative difference of 0.033±0.063 between grc4f and PYTHIA is found. In Table ^ 
we conservatively take the statistical error on this difference as systematic uncertainty for both 
the “eey*” and the “eeZ” region. 

Uncertainties in the description of the background processes cause systematic errors on the 
background subtraction. The contributions to these errors coming from e’''e“ —> qq, two-photon 
processes and other four-fermion backgrounds have been tested separately for each background. 
By inverting or removing one or two selection cuts while the others were left unchanged, we pro¬ 
duced samples where the analysed background process is signihcantly dominant. We obtained 
the systematic uncertainty for the various background processes by taking the relative difference 
of the total event numbers in data and Monte Carlo after the statistical error on the data has 
been subtracted quadratically. The relative systematic uncertainties we obtain are 0.17 for the 
e“''e“ —>■ qq background, 0.14 for the four-fermion backgrounds and 0.24 for the two-photon pro¬ 
cess e'''e“ —> e’''e“-|-hadrons. The numbers of background events after all cuts are 1.8±0.5±0.3, 
1.4±0.1±0.2 and 3.1±0.8±0.7 for e'''e“ —qq, e’''e“ —>• 4f and e'''e“ —>• e'''e“ -f hadrons, respec¬ 
tively. Since other background sources are much smaller, their systematic uncertainties have 
been neglected. The total error on the background events corresponds to a relative systematic 
uncertainty for the obtained cross-sections of 0.070 for “eej*” and 0.018 for “eeZ” (see Table 
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8 Results and conclusion 


Using a data set with an integrated luminosity of (54.7±0.1±0.2) pb~^ collected with the 
OPAL detector at y/s = 183 GeV, 27 candidate events for the process e’''e“ —> (e)eZ/ 7 * have 
been selected, with an expected background of 6.5±0.9±0.8 events, where the first errors are 
statistical and the second are systematic. The expected number of signal events is 16.7±0.3±1.9. 


From the observed number of events, the cross-section for this process has been calculated 
separately for a high mass, “eeZ”-like, and a low mass, “ee 7 *”-hke region, with a cut at a 
hadronic mass of 60 GeV. Within the phase space of our kinematic signal dehnition described 
in Section we measure a cross-section of (4.1±1.6±0.6) pb for “ee 7 *”-like events, and a cross- 
section of (0.9±0.3±0.1) pb for “eeZ”-like events, including the systematic errors obtained in 
Section ^ 


To investigate in more detail the kinematic properties of the Gompton scattering process, 
we show in Fig. ^ the differential distributions of the Mandelstam variables from Equation [^, 
and of the cosine of the scattering angle 6* of the Z/ 7 * in the ey rest system with respect 
to the incoming 7 direction. All quantities have been obtained from the measured momenta 
and energies of the electron and the hadronic system, and are compared with the grc4f signal 
prediction. The structure observed in the \/§ distribution is consistent with expectation for the 
7 * (low ^/§) and Z (high a/J) contributions. An excited electron, contributing in the s-channel 
of the signal diagram, would show up as a peak at the e* mass. For illustration, we have 
added in Fig. as a dashed histogram the expected contribution of events with an excited 
electron of mass m* = 120 GeV, normalized to a product of cross-section times branching ratio 
of 5 pb. It has been obtained by applying our analysis to fully simulated events of the process 

m. 


e^e 


ee* —*• eeZ from the Monte Garlo generator EXOTIG 


In the t and u distributions, peaks near zero are predicted, with a long tail to large negative 
values in the case of t. The measured t distribution agrees well with this prediction and 
also shows the drop near zero, due to angular acceptance cuts. The observed u distribution 
seems somewhat broader than predicted, though it is still consistent within errors. The entries 
at positive values are due to detector resolution. The distribution of cos 6 ** peaks at —1, 
corresponding to small u which is typical for Gompton processes. 


In summary, we have reported the first observation of Z/ 7 * production in Gompton scat¬ 
tering of quasi-real photons which is a subprocess of the reaction e’''e“ —(e)eZ/ 7 *. The pre¬ 
dictions of the grc4f and PYTHIA Monte Garlo programs, as listed in Table are both in 
good agreement with our results. The data statistics are not yet sufficient to be sensitive to 
the difference between the predictions of the Monte Garlo generators. 
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Figure 2: Distributions of m^q and after the preselection. The histograms show the contri¬ 
butions from the various Monte Carlo simulations and the points are the data. The sharp lower 
edges of the spectra are due to preselection cuts. Only statistical errors are shown. 
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Figure 3: Distributions of the selection variables. The histograms show the various Monte Carlo 
simulations and the points are the data. The arrows point into the selected regions. The cuts 
have been applied successively. Only statistical errors are shown. 
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data 


background 



E, (GeV) (GeV) 

Figure 4: Distributions of niqq and Eg after all cuts. In Figures (a) and (c) the signal is 
simulated using the grcff generator. In (b) and (d) we have used PYTHIA instead of grcff. 
The histograms show the various Monte Carlo simulations and the points are the data. Only 
statistical errors are shown. 
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OPAL 


data 


background 



Figure 5: Distributions of t, u and cos 9* after all cuts. The histograms show the 

various Monte Carlo simulations and the points are the data. The dashed histogram in (a) is a 
hypothetical e* signal, as described in the text. Only statistical errors are shown. 
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